Astronomy & Astrophysics manuscript no. 10187 
July 21, 2008 



© ESO 2008 



Amplitude variability and multiple frequencies in 44 Tau: 

2000 - 2006 

M. Breger and P. Lenz 



00 

o 
o 

(N 



6 



Institut fiir Astronomie der Universitat Wien, Turkenschanzstr. 17, A-1180 Wien, Austria 
e-mail: michel.breger@univie.ac.at 



Received date; accepted date 



ABSTRACT 



rn 

O 
oo 
O 



X 



Aims. This study has three principal aims: (i) to increase the number of detected pulsation modes of 44 Tau, especially outside the 
previously known frequency ranges, (ii) to study the amplitude variability and its systematics, and (iii) to examine the combination 
frequencies. 

Methods. During the 2004/5 and 2005/6 observing seasons, high-precision photometry was obtained with the Vienna Automatic 
Photoelectric Telescope in Arizona during 52 nights. Together with previous campaigns, a data base from 2000 to 2006 was available 
for multifrequency analyses. 

Results. Forty-nine pulsation frequencies are detected, of which 15 are independent pulsation modes and 34 combination frequencies 
or harmonics. The newly found gravity mode at 5.30 cd"' extends the known frequency range of instability. Strong amplitude variabil- 
ity from year to year is found for the t = I modes, while the two radial modes have essentially constant amplitudes. Possible origins 
of the amplitude variability of the £ = I modes, such as precession of the pulsation axis, beating and resonance effects are considered. 
The amplitudes of the combination frequencies, / -I- fj, mirror the variations in the parent modes. The combination parameter, which 
relates the amplitudes of the combination frequencies to those of the parent modes, is found to be different for different parents. 
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Q\^ 1. Introduction 



The asteroseismic comparison between measurements and the- 
oretical models of the radial and nonradial pulsation of stars on 
and near the main sequence involves the exact matching of the 
values of the pulsation frequencies, as well as reproducing the 
frequency range in which instability occurs. However, it is ob- 
servationally difficult to determine the important borders of the 
frequency range in which a particular star is unstable. Due to 
the small amplitude growth rates of modes on the edges of the 
frequency instability range, only modes with small photometric 
amplitudes are expected to appear in these frequency regions. 

If we consider the 6 Scuti star FG Vir as an example (Breger 
et al. 2005), all photometric amplitudes in the 30 - 45 cd"' region 
had amplitudes less than 1 .0 millimag. In fact, only the addition 
of 2 1 8 additional nights of photometry to the previous data al- 
lowed the detection of the pulsation modes near 40 cd '. In ad- 
dition to the problem of small amplitudes in specific frequency 
regions, it is also essential to separate the 'real' modes from 
the combination frequencies, /J + fj, which occur in these high- 
frequency and low-frequency regions. This separation requires 
very high frequency resolution so that the chance of accidental 
agreements between the expected and measured frequency val- 
ues is essentially reduced to zero. In practice this requires cam- 
paigns covering several years. 

The requirement of multiyear campaigns is strengthened 
by the amplitudes of many pulsators showing variability on 
timescales of months and, especially, years. To determine what 
types of modes show amplitude variability and to search for 
the astrophysical origin of this amplitude (and frequency) vari- 
ability requires extensive observing campaigns. Finally, success- 
ful asteroseismology also demands the identification of individ- 



ual pulsation modes. Both spectroscopic and photometric tech- 
niques have been applied: however, the large amounts of observ- 
ing time required are easier to obtain on smaller telescopes with 
photometers. 

The (5 Scuti star 44 Tau is ideal for a detailed study because 
of its extremely low rotational velocity of 3 + 2 km s"' (Zima 
et al. 2007), for which second-order effects of rotation can be 
neglected, and because its amplitudes place the star between the 
high-amplitude pulsators (HADS) with dominant radial modes 
and the average low-amplitude 6 Scuti variable with mostly non- 
radial modes. Zima et al. (2007) present high-resolution spec- 
troscopy of 44 Tau, which was used together with radial velocity 
data to derive the photospheric element abundances and provide 
m values for the pulsation mode identifications. Six axisymmet- 
ric and two prograde modes were identified. 

In their asteroseismic study of 44 Tau, Lenz et al. (2008) 
present spherical mode degree identifications for the dominant 
modes of 44 Tau. Photometric amplitude ratios and phase differ- 
ences suggest two radial, four £ - I, and three { - 2 modes. Due 
to the measured \ogg value of 3.6 + 0.1, both main sequence 
and post-main sequence models were examined. The identified 
modes can be fitted well in both evolutionary stages. The pre- 
dicted frequency ranges for unstable modes are in good agree- 
ment with the observed ranges adopted by Lenz et al. (2008). 
Due to the different envelope structure of the main sequence 
and post-main sequence model, the borders of instability differ 
slightly. Consequently, new independent frequencies in the low- 
frequency and high-frequency regions possess important infor- 
mation for improving the models. 

We have previously presented 90 nights of photomet- 
ric observations obtained during the three observing seasons 
2000/2001, 2001/2002, and 2002/3 (Antoci et al. 2007: this pa- 
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per also lists a detailed observational history of the star). The 
data from three years led to detecting 29 frequencies, of which 
13 were independent pulsation modes. The two comparison stars 
used showed low-frequency variability at the millimag level, 
which is, regrettably, quite common. Otherwise, these measure- 
ments are very accurate. Consequently, the publication omitted 
the to 5 cd"' region awaiting future measurements with differ- 
ent comparison stars, as reported in the present paper. 



2. New photometry obtained 2004/5 and 2005/6 

The APT measurements were obtained with the T6 0.75 m 
Vienna Automatic Photoelectric Telescope (APT), situated at 
Washington Camp in Arizona (Strassmeier et al. 1997, Breger & 
Hiesberger 1999). The telescope has been used before for several 
lengthy campaigns of the Delta Scuti Network, which confirmed 
the longterm stability and millimag precision of the APT pho- 
tometry. The V and y filters of the Stromgren system were used. 
Only high-quality data were retained with the following number 
of nights and hours: 2004/5 observing season - 143 hours of pho- 
tometry during 26 nights, 2005/6 observing season - 173 hours 
also during 26 nights. 

The rapid position changes of the APT made it possible to 
observe 4 comparison stars: HD 23626 (called CI, spectral type 
GO), HD 24301 (called C2, GOIV), and the two stars used dur- 
ing previous campaigns HD 25867 (FIV) and HD 25768 (F8). 
The new comparison stars confirmed the suspicion that the two 
previously used comparison stars are slightly variable with pe- 
riods longer than 1 day, which would be typical of Gamma 
Doradus stars found at these spectral types. The results were 
already demonstrated in Fig. 3 of Breger (2007) as a warning 
to observers about common millimag variables. These results 
also presented evidence for 0.936 and 0.97 cd ' frequencies in 
HD 25867 and 0.885 cd' in HD 25768 with millimag ampli- 
tudes. The results allowed us to rereduce the 2000/1, 2001/2, 
and 2002/3 data by prewhitening these frequencies from the 
comparison-star data. No variability was detected for the two 
new comparison stars and the 2004/5 and 2005/6 measurements 
were reduced by interpolating between these comparison stars. 
The comparison stars agreed to + 2 mmag per single measure- 
ment in both the y and v passbands. 

The new light curves with the final fits are not shown, since 
the diagrams look similar to the 2000 - 2003 measurements of 
44 Tau published by Antoci et al. (2007). The agreement be- 
tween the observations and fits are excellent with residuals per 
single measurement of + 2 mmag in the >■ and + 3 mmag in the 
V passbands. 



3. Frequency analyses 

The pulsation frequency analyses were performed with a pack- 
age of computer programs with single-frequency and multiple- 
frequency techniques (PERIOD04, Lenz & Breger 2005), which 
utilize both Fourier and multiple-least-squares algorithms. The 
latter technique fits up to several hundreds of simultaneous si- 
nusoidal variations in the magnitude domain and does not rely 
on sequential prewhitening. The amplitudes and phases of all 
modes/frequencies are determined by minimizing the residuals 
between the measurements and the fit. The frequencies can also 
be improved at the same time. 

New frequencies of pulsation were searched for by applying 
Fourier analyses to the data prewhitened by the previously de- 
tected frequencies. However, due to strong amplitude variability 
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Fig. 1. Fourier spectra of the residuals in the 5-35 cd"' region 
after prewhitening 49 frequencies. The drawn curve represents 
the calculated significance level. 



from year to year, as well as the different amplitudes and phases 
found in the v and y passbands, for each year separate multifre- 
quency solutions were made, which could then be prewhitened 
to search for additional frequencies. This approach works well 
for modes with amplitudes in excess of 1 mmag. For the many 
modes with small amplitudes, we improved the signal-to-noise 
ratios by combining the data from the v and y passbands. This 
was done by scaling the v data by an experimentally determined 
factor of 0.69 and increasing the weights of the scaled data ac- 
cordingly. In principle, this could create problems due to small 
phase shifts, but this is negligible for modes with such small am- 
plitudes. Following Breger et al. (1993), a significance criterion 
of an amplitude signal-to-noise ratio of 4.0 (which corresponds 
to a power signal-to-noise ratio ~12.6) was adopted for all fre- 
quencies. 

The results are shown in Table 1. For /12 to fug, we only list 
the combined amplitude for the v and y passbands. Due to the 
scaling of the v data, this corresponds to the y passband, whose 
effective wavelength is similar to that of the standard Johnson 
V passband. The pulsation mode identifications listed are from 
Zima et al. (2007) and Lenz et al. (2008). The formal uncertain- 
ties in the amplitudes are listed in Table 2. These values need to 
be considered in discussing the amplitude variability. 

3.1. The main pulsation region 

The main pulsation region, in which we find independent modes, 
ranges from 5.3 to 12.7 cd"', which is more than the previously 
known range. Two previously unknown frequencies, represent- 
ing independent pulsation modes, were discovered. These are 
found at 5.30 and 7.79 cd"', respectively. In both cases, the am- 
plitude signal-to-noise ratio is higher than 10, representing cer- 
tain detections. The 7.79 cd"' mode is near a 1 cd"' alias of the 
previously known 6.80 cd"' mode. However, we found it impos- 
sible to identify the new frequency as an artifact of the 6.80 cd"' 
mode, caused by amplitude variability. 

Are there additional modes seen in the data? Two more peaks 
are near the statistical detectability limit of the amplitude signal- 
to-noise ratio of 4.0: a peak at 21.243 cd"' has an amplitude 
near 0.15 mmag (Fig. 1), while the 3.92 cd"' peak is seen only 
in the 2004/5 and 2005/6 data with a height of 0.4 mmag. The 
difficulty does not lie with the signal, but with the noise figure. 
We do not regard the calculated noise as reliable; in these fre- 
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Table 1. Frequencies and amplitudes in the Stromgren v and y filters. 





Frequency 


Mode 








Amplitude 


in minima 


g 










cd-' 


{t,m) 


2000/1 


2001/2 


2002/3 


2004/5 


2005/6 






or ID 


V 


y 


V 


y 


V 


y 


V 


y 


V 


y 


/l 


6.8980 


(0,0) 


59.64 


27.37 


39.41 


27.46 


39.32 


26.94 


39.23 


27.42 


39.53 


27.40 


/2 


7.0060 


(1,1) 


19.11 


13.30 


16.80 


11.86 


13.86 


9.19 


10.00 


6.92 


8.29 


5.64 


/3 


9.1174 


(1,1) 


16.76 


11.48 


21.09 


14.56 


17.62 


11.91 


4.87 


3.37 


0.81 


0.61 


/4 


11.5196 


(1,0) 


18.16 


12.73 


16.53 


11.65 


16.72 


11.70 


10.25 


7.21 


8.78 


6.36 


/5 


8.9606 


(0,0) 


13.81 


9.53 


13.75 


9.24 
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0.15 
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0.32 
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Note: Some close frequency pairs are near the annual frequency resolution. Consequently, annual /31 values are not given and for 2002/3, some 
amplitudes marked with a colon are uncertain. 



quency ranges, there are many correctly identified combination 
frequencies, which have been prewhitened. However, after each 
prewhitening, a zero amplitude remains at the prewhitened fre- 
quency value, while there should actually be some noise with 



different phasing left. After prewhitening a number of these fre- 
quencies, the overall noise is reduced. Evidence of this effect 
can be seen in Fig. 2, since the expected increase in noise to- 
ward lower frequencies is not seen. Consequently, modes near 
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2004/5 - 2005/6 







1 2 3 
Frequency (c/d) 


4 5 


0.4 


- 


1 1 1 1 1 
2000 to 2006, after prewhitening 




0.2 
0.0 


i 


d^^l^flil^lt^ 



Fig. 2. Fourier spectra of the low-frequency region of 44 Tau. 
Note the many identified combination frequencies and their 1 
cd ' aliases. We have shown the Fourier spectrum of the final 
residuals for both the 2004/5 and 2005/6, which used excellent 
comparison stars, and the total data set. The statistical signifi- 
cance levels are also shown. 
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Fig. 3. Two combination frequencies with similar values. The 
diagram shows that the frequency resolution of the six-year data 
set is high enough to resolve these two frequencies. 



the formal statistical detection limit should be treated with cau- 
tion; i.e., we have not accepted the two peaks near the detection 
limit. An additional argument can be made that such an uncer- 
tain detection should be seen in two or more independent data 
sets: this is not the case here. The two peaks may, therefore, be 
caused by noise. 

What mode does the new frequency at 5.30 cd ' represent? 
The mode at 5.30 cd ' increases the range of excited frequencies 
into the g-mode domain, since the radial fundamental mode is 
found at 6.90 cd ' . The identification of the 6.90 cd ' as the ra- 
dial fundamental mode is secure because of photometric, as well 
as spectroscopic, mode identifications, the observed frequency 
ratio of 0.770 between the observed 6.90 and 8.96 cd"' modes 
and the known parallax. The frequency of the newly discovered 
5.30 cd"' mode also gives a frequency ratio of 0.769 relative to 
the known radial fundamental mode: the near agreement with 
the radial ratio of 0.770 serves as a warning that the frequency 
ratio alone is not a sufficient criterion for mode identification. 

3.2. Combination frequencies 

The excellent frequency resolution provided by the coverage 
from 2000-2006 makes it possible to accurately identify those 
peaks that are combinations, fi+fj, as well as 2/i harmonics. We 
find that the peaks outside the 5.3 to 12.7 cd"' range can be iden- 
tified as combinations or harmonics. We find it remarkable that 
even some triple combinations are seen. The parent modes with 
the highest amplitudes were involved most often in the combi- 
nation frequencies, which is not surprising. 

An interesting situation arises because of the similarity in 
some predicted frequency combinations. Numerically, the fre- 
quency values of three { = I modes /2, /4, /e, and the radial 
mode fi are related in the sense that 

/2 + /4 =/5+/6 + 0.0039 cd-'. 
This is a very close, probably accidental, agreement. The small 
frequency separation is only of the order of the aliasing caused 
by the gaps between the annual observing seasons. For the whole 
time period of six years, the frequency resolution is high enough 
to separate both combinations. Indeed, Fig. 3 shows that both 
are present, but not at equal amplitudes. However, for a single 
observing season, the frequency resolution is too low to separate 
the two combination frequencies; consequently, we only derive 
a combined amplitude. 

A similar situation exists for combination frequencies 
formed by frequency differences /s - /a vs. /4 - /g near 1.95 cd"' 
and /e - /a and /4 - /s near 2.56 cd"'. However, in these cases 
the power spectrum covering five years only shows unique, sin- 
gle detections. 



Table 2. Formal uncertainties of the amplitude values in mmag 



Observing season 


o-(v) 


o-(y) 


cr(combined >■) 


2000/1 


0.07 


0.06 


0.04 


2001/2 


0.12 


0.12 


0.07 


2002/3 


0.19 


0.15 


0.10 


2004/5 


0.09 


0.07 


0.05 


2005/6 


0.08 


0.07 


0.05 



Extensive Monte Carlo simulations give larger uncertainties 
with a factor of ~ 1.35. 
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4. Amplitude variability: Observations 

It has already been noticed by Poretti et al. (1992) that in 44 
Tau the amplitudes of /2 and /g are variable. The new exten- 
sive data covering five observing seasons in six years allow us 
to examine these and other variations in more detail. For the 
dominant modes with mode identifications, the amplitude vari- 
ability is shown in Fig. 4. Here the long 2000/1 observing sea- 
son was subdivided into two parts with the average time strings 
spaced 40 days apart. Separate multifrequency solutions show 
that the small variations are as expected from the larger year- 
to-year variations and rapid changes are not seen. Note that the 
steady amplitude increase in f(, after 2004/5 is confirmed when 
the 2005/6 season is subdivided into two parts. 

The observed amplitude variability can be summarized as 
follows: 

(i) The two radial modes at /i and fs show very little or 
no amplitude variability from 2000 to 2006. For both modes, 
the new amplitude values also agree with the amplitudes from 
1989/1990 given by Poretti et al. (1992). The results are consis- 
tent with amplitude constancy. 

(ii) All four i - \ modes (/2, /s, /4, and /s) show strong 
amplitude variability with decreasing amplitudes between 2000 
and 2006. The timescale of variation is several years with little 
variation seen within an observing season. A general trend of 
decreasing amplitudes from 2000 to 2006 exists. However, the 
different modes do not vary together: the differences in behavior 
of the four modes considerably exceed the observational uncer- 
tainties. The different behavior is not connected to the fact that 
the modes at /2 and /3 are prograde, while fn is axisymmetric. 

(iii) The { -2 modes have smaller photometric amplitudes, 
but the amplitude variability exceeds the statistical uncertainties. 
Otherwise, no systematic behavior is seen. 

4.1. Amplitude variability of the combination frequencies 

The observed combination frequencies in 6 Scuti stars may be 
created by two mechanisms. Resonant excitation of a pulsation 
mode close to the combination frequency leads to frequency syn- 
chronization at exactly fi + fj. The second possibility is that a 
nonlinear response of the medium in the outer layers of the star 
causes pulsation signals at combination frequencies. The latter 
mechanism is likely to be responsible for combination frequen- 
cies in DA and DB white dwarfs (see e.g. van Kerkwijk et al. 
2000). 

If the same mechanism is responsible for the amplitude vari- 
ability in 44 Tau, we can expect the combination frequencies to 
mirror the amplitude variations of the parent modes. The strong 
amplitude variability of the { = I modes makes it possible to 
examine the amplitudes of the combination frequencies, fi + fj, 
involving one or two of these modes. Three difficulties arise: if 
both fi and fj are variable, it becomes difficult to separate the 
effects from the two components. This can be avoided if one of 
the parents of the combination frequencies has a constant ampli- 
tude; i.e., in 44 Tau one of the parents should be a radial mode. 
Secondly, the parent modes need to have high amplitudes for 
a combination frequency to be visible. This restricts us to the 
dominant modes. Thirdly, the amplitudes of the combination fre- 
quencies are small and affected by random or systematic noise; 
consequently, we have omitted the low-frequency region. The 
combination frequencies shown in Fig. 5 obey the criteria and a 
very good correlation is seen. 

The agreement between the amplitudes of the combination 
and the parent modes is very good. This argues against an origin 



of the combination frequencies as independent modes excited 
by resonance: with resonance the power can be redistributed 
between the three modes, but not just decreased. The measure- 
ments, however, match the effects due to nonlinear mixing of the 
eigenmodes, as discussed by Garrido & Rodriguez (1996). 

The amplitudes given in Table 1 are used to calculate the 
combination parameter, //, which relates the amplitude of the 
combination frequency with those of the parents. 



"^comb 



■ fi-Ar Aj 



(1) 



where A, and Aj are the amplitudes of the parent modes. This 
definition is similar to the definition of van Kerkwijk et al. 
(2000). We have calculated the sizes of the combination param- 
eter for the different combinations including a complete treat- 
ment of the propagation errors. The uncertainties in the ampli- 
tudes (given in Table 2) were used after multiplying the values 
by a factor of 1.35 derived from Monte Carlo simulations us- 
ing the real data. The data from the different years were aver- 
aged with weights derived from the known uncertainties of the 
annual amplitudes. The resulting coefficients, fi, are shown in 
Table 3. Regrettably, the coefficients of other interesting combi- 
nations are too uncertain to be listed. 

Table 3. Relationship between the amplitudes of combination 
frequencies and parent modes. 



Combination frequencies 

fi + fj {t,m) + {{,m) 



Combination parameter 

y" = A,„„,b/(A,- ■ Aj) 



A+fs 

/1+/2 
/1+/3 
/1+/4 

/5+/4 

Uncertain values 

fi+fe 

/5+/2 



(0,0) + (0,0) 

(0,0) + (1,1) 
(0,0) + (1,1) 
(0,0) + (1,0) 
(0,0) + (1,0) 

(1,0) + (1,?) 
(0,0) + (1,1) 



0.0027 ± 0.0005 

0.0028 ± 0.0004 
0.0028 ± 0.0002 
0.0034 ± 0.0003 
0.0092 ± 0.0005 



0.0045 : 
0.0071 : 



0.0018 
0.0019 



One may assume that the combination frequency involving 
a radial and a ^ = 1 mode shows the same geometrical cancel- 
lation effects as the { = 1 mode of the combination, while there 
is no light cancellation for the radial mode. Consequently, the 
combination parameters in Table 3 can be considered to already 
be free of geometrical cancellation effects. 

We find that 

(i) the combinations involving the radial fundamental mode, 
/i, all have similar combination parameters near 0.003. This in- 
cludes the combinations with the { - I and { -0 modes. 

(ii) The two radial modes have different combination param- 
eters; although the first overtone, fi, has a considerably smaller 
amplitude than the fundamental mode, /i , the combination fre- 
quencies (/1-1-/4) and ifs+fA) have similar amplitudes, Acomb- 

(iii) Some (weak) evidence exists that the combinations of 
the first radial overtone, fs, with other nonradial modes also have 
higher values of the combination parameter. 

5. Theoretical interpretation of the observed 
amplitude variability of the i =1 modes 

The exceptionally slow rotation of 44 Tau needs to be consid- 
ered when analyzing amplitude variability. The reason for the 
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Fig. 4. Amplitude variability of radial and nonradial modes found in 44 Tau for the time period 2000-2006. The amplitudes in both 
the y (~ V) and the v passbands were used after scaling the v amplitudes. The eiTor bars are shown and are usually the same size as 
the symbols used. Note the relative constancy of the radial modes and the large variability of the ^ = 1 modes. The variations within 
the 2000/1 and 2005/6 (/g only) observing seasons confirm the long-term trends. 



slow rotation of 44 Tau has not yet been revealed. Common rea- 
sons for slow rotation are magnetic braking, a close companion 
star, or planets. However, the observed spectra do not show any 
indications of a strong magnetic field or a binary star. 

We examined different scenarios to explain the observed am- 
plitude variability in 44 Tau. Only for one of the four /" = 1 
modes, /6, was a full modulation period presumably covered. Its 
duration was approximately 5 years. We now consider several 
mechanisms that cause amplitude modulation: beating of close 
frequencies, resonant coupling of modes, and/or the precession 
of the pulsation axis that causes a variation in the visibility of a 
pulsation mode. 

5.1. Amplitude variability due to precession of the pulsation 
axis 

A possible scenario that causes amplitude variability for { - I 
modes, but does not influence the photometric amplitudes of ra- 
dial modes, is the precession of the pulsation axis, which is usu- 
ally assumed to coincide with the rotation axis. The variation 
in the stellar inclination angle results in a changing cancellation 
of light for { - \ modes, while the visibiUty of ^ = modes 
remains the same. 

Zima et al. (2007) have identified two ( - \ modes, /2 
and /3, as prograde modes {m - 1) and /4 as an axisym- 
metric mode (m - 0). For the fourth £ - \ mode, /e, the 
azimuthal number could not be identified. The measured in- 
clination in 2004 was 60 + 25°. We computed mode ampli- 
tudes in the Stromgren y band according to the equations given 
in Daszyhska-Daszkiewicz et al. (2003) assuming the intrinsic 
mode amplitude, e, to be 0.000125. This value was estimated 



by scaling the predicted photometric amplitudes to the observed 
amplitudes. The actual value for s is unimportant in this study 
because we only want to examine the effect of the inclination 
on mode visibility and not to predict absolute photometric am- 
plitudes. Our computation also takes limb darkening effects into 
account. The frequency dependence between the relative grav- 
ity changes and the radius changes needs to be considered, and 
one has to consider the frequency of the given mode as well. 
However, the influence of the frequency value on the visibility is 
weak. In this examination we relied on the post-main sequence 
model presented in Lenz et al. (2008). It should be noted that the 
specific model parameters do not change the main results in this 
section significantly. 

The dependence of photometric amplitudes on stellar in- 
clination in the case of the ( - \ mode /g is shown in Fig. 
6. It can be clearly seen that the visibility of an axisymmet- 
ric and non-axisymmetric mode is different. The intrinsic mode 
amplitudes are assumed to be constant, and the different visi- 
bility is only due to cancellation effects. The amplitude of a 
((,m) - ( 1 , 0) mode decreases with increasing inclination and 
finally is fully cancelled out at 90° (equator-on view). The vis- 
ibility of a {{,m) = (1, 1) mode is, however, highest for large 
inclinations and lowest at small inclinations (pole-on view). 

The effect of the precession of the pulsation axis on the vis- 
ibility of the { - \ modes /2 and /& is shown in Fig. |2l Here 
we assume that the inclination varies between 30 and 90° and 
that /e is an axisymmetric mode. This diagram looks very simi- 
lar when /3 and f^ are plotted. The predicted change of visibil- 
ity of the axisymmetric /g mode mimics the observed amplitude 
variability very well. This would indicate an actual inclination 
angle close to 85° for the year 2004. The inclination angle was 
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Fig. 5. The relation between amplitude variability of combination frequencies and their parents. We have selected only modes for 
which one of the parents is a radial mode with constant (or nearly constant) amplitudes. The diagram demonstrates the excellent 
correlations. 
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Fig. 6. Dependence of mode visibility in the Stromgren y filter 
on the stellar inclination shown for an axisymmetric and a non- 
axisymmetric f = 1 mode with the same frequency. The shaded 
area marks the measured limits for the inclination of 44 Tau. 
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Fig. 7. Predicted amplitude variability due to precession of the 
pulsation symmetry axis for an inclination varying between 30 
and 90°. Lower panel: assumed variation of the inclination in an 
arbitrary time unit. Upper panel: Stromgren y amplitudes of the 
(i,m) = ( 1 , 0) mode /g and the (^, m) = ( 1 , 1 ) mode /i computed 
for the corresponding inclination values. 



previously determined in 2004 with a mean value of 60 + 25°. 
Therefore, an inclination close to 85° is still within the given un- 
certainty limit. However, the behavior of the non-axisymmetric 
mode /2 is clearly different and its amplitude expected to be 
highest at large inclination angles. As can be seen from Fig. [H 
this is not observed. Furthermore, the second axisymmetric { = I 
mode, /4, does not go to zero amplitude simultaneously with f^. 



Consequently, the precession of the stellar inclination cannot ex- 
plain the observed amplitude variability of all the { - I modes. 
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5.2. Beating effects 

Beating effects between close frequencies often explain short- 
scale amplitude variability in 6 Scuti stars (e.g., FG Vir, Breger 
& Pamyatnykh 2006). Since 44 Tau rotates at a very low rate 
(3 + 2 km s"'), the rotationally split components are separated 
by only a very small amount. In such a case the components of 
the multiplet may cause amplitude modulation if the frequencies 
are too close to be resolved in the observed data set. 

Since we intend to examine this possibility for all observed 
{ = 1 modes in 44 Tau, we carried out tests with artificial data 
sets to determine the expected amplitude and phase modulation 
for three-mode beating. The differing visibility of the compo- 
nents of the multiplet depending on the inclination angle needs 
to be taken into account. Consequently, we examined the situ- 
ation that occurs for the two boundary cases for the inclination 
of 44 Tau, / = 35 and 85°. The inset diagram in Fig. |8] shows 
the visibility of the multiplet components for these inclination 
angles. For an inclination of 85°, the observed amplitude of the 
central peak (m = 0) is close to zero, and we essentially have a 
similar situation as for two-mode beating. 
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Fig. 8. Expected amplitude and phase modulation for three- 
mode beating considering geometrical visibility conditions for 
an inclination angle of 35° (left panels) and 85° (right panels). 
The amplitudes of the components of the triplet are given in the 
inset diagrams. The rotational splitting was computed for the 
measured mean rotation rate of 3 km s"' . 



To examine the effect of three-mode beating, we computed 
the rotational splitting of the t - I mode /(, according to 2nd- 
order theory assuming a rotation rate of 3 km s"' . The frequency 
separation of 0.01 1 cd"' is almost equidistant at such a low rota- 
tion. We created an artificial time string using the corresponding 
amplitudes for the given limits of the inclination angle. This data 
set was then analyzed with PERIOD04 assuming a single unre- 
solved frequency centered at the axisymmetric mode. The am- 
plitude and phase variability that occurs in such a case is given 
in Fig.[8]for the inclination angles 35 and 85°, respectively. 

An important observable is the occurrence of a phase change 
at minimum amplitude. The minimum of the amplitude modu- 
lation cycle has been observed only for /e. No significant phase 
change was detected. However, the determination is uncertain 
because the exact time of the minimum is not known, and due 
to observing gaps, the uncertainties in the phase determination 



are large. Consequently, we can only draw conclusions from the 
expected timescales of the beating effect. 

The period of amplitude modulation, Pnmd, corresponds to 
l/Avo,±i where Avo-ti is the mean frequency spacing between 
the m = and \m\ = 1 components. As already mentioned ear- 
lier, at an inclination angle of 90°, the observable amplitude of 
the axisymmetric mode vanishes and two-mode beating occurs 
with a modulation period 1/Av_i +i, where Av_i +i denotes the 
frequency separation between the two non-axisymmetric com- 
ponents of the triplet. Since the inclination of 44 Tau may be as 
high as 85°, we cannot exclude such a case. 

These results allow us to examine the expected timescales 
of beating effects for 44 Tau. At the upper limit of the rotation 
rate, 5 km s"', the rotational splitting amounts to 0.0185 cd"', 
which causes amplitude modulation periods of 54 d for three- 
mode beating and 27 d for two-mode beating. At the lower limit 
of the rotation rate, 1 km s"' , the rotational splitting of a trapped 
^ = 1 triplet is 0.00383 cd"'. This corresponds to a beat period 
of 260 d for three-mode beating and 130 d for two-mode beat- 
ing between the non-axisymmetric modes at the high inclination 
limit. It can be clearly seen that the timescales for amplitude 
modulation due to beating between the components of a ^ = 1 
multiplet are too short within the measured limits of rotation. 
Beating timescales, that would match the observed variability, 
require an improbable rotational velocity of <0.1 km s '. 

We conclude that beating effects between components of a 
rotationally split triplet cannot explain the observed behavior. In 
principle, the observed long timescales could also be caused by 
a very close mode with a different spherical degree. The obser- 
vational test reUes critically on the observed phase shift at the 
time of minimum amplitude. Unfortunately, there is not enough 
data to carry out this test. 



5.3. Resonance effects 

Another cause for amplitude modulation is the resonant coupling 
of pulsation modes (e.g., Nowakowski 2005). Unstable pulsa- 
tion modes are likely to be coupled with linearly stable modes 
of high spherical degrees. These modes may not be detected by 
photometric and spectroscopic methods due to cancellation ef- 
fects across the stellar disk. However, to examine whether the 
amplitude variability of the { = I modes is due to resonance, the 
expected timescales of amplitude modulation can be compared 
to the observed variation. 

Resonance effects show periodicities on timescales of the in- 
verse linear amplitude growth rate of the unstable mode (e.g., 
Moskalik 1985). Table |4] lists the expected timescales, Tres, for 
all four observed £ = I modes, according to a pulsation model 
for 44 Tau. 



Table 4. Predicted timescales for resonance effects. 



Frequency [cd '] 


rres ~ 1/r [yrs] 




7-1:1 


1 [yrs] 








V,-„ 


= lkms-' 


* rot 


=5 km s-i 


/> 7.01 


122.4 




137 




4.4 


/3 9.12 


90.5 




55 




2.6 


fe 9.56 


16.1 




>150 




17.1 


/4 11.52 


1.3 




45.6 




1.8 



M. Breger and P. Lenz: Amplitude variability and multiple frequencies in 44 Tau: 2000 - 2006 



The observed variability of the amphtude of the four £ - I 
modes occurs on a similar timescale, while the predicted modu- 
lation periods, Ties, are very different. In particular, the observed 
variation in the amplitudes of /2 and fj, has a significantly shorter 
timescale than predicted. On the other hand, for /4 the theoretical 
modulation timescale is too short. 

We also investigated the possibility of a 1:1:1 resonance be- 
tween the frequencies of an /" = 1 multiplet (see e.g. Buchler, 
Goupil & Hansen 1997). The modulation period is riii w 
|l/(v+i +V-1 -2vo)| where the subscript denotes the m value. For 
a model of 44 Tau with a rotational velocity at the lower limit 
of the measured rotation rate, the predicted values for tiii have 
reasonable timescales with the exception of /& (see Table|4]l. For 
the upper limit of the rotation rate of 5.0 km s"', the computed 
timescales are too short for /2, fj, and /4 and still too long to 
explain the observed variability of /e. 

Another difficulty concerns the nondetection of the other 
components of the rotational triplets. The components of the 
triplet are separated by ~ 0.02 cd ' at the upper limit of the stel- 
lar rotation rate. We have only detected a single possible case 
(9.58 cd"') with a small amplitude: regrettably, this mode does 
not have a mode identification. The missing detection of multi- 
plet structures provides an argument against the 1:1:1 resonance 
hypothesis for rotation rates near our upper limit: due to the ob- 
served decrease in the amplitude of the observed mode, the am- 
plitudes of the other components of the multiplet should already 
be large enough to be detectable. Consequently, while resonance 
between the components of a multiplet may explain the ampli- 
tude variability of a fraction of the { = I modes, the 5-yr am- 
plitude modulation of /g cannot be explained by the 1:1:1 reso- 
nance. 

This section has shown that resonance may be responsi- 
ble for the observed amplitude variability of specific modes. 
However, for definitive conclusions a longer time baseline for 
the observations would be required. Such data would also clar- 
ify whether the fascinating observation of a near-simultaneous 
decrease of the amplitudes of all the observed £ = I modes is 
accidental. 



6. Impact of the new frequencies on the pulsation 
models of 44 Tau 

One of the main goals in this study was to search for new sig- 
nificant frequencies that extend the previously known frequency 
range of 6.34 - 12.70 cd ' and to compare the new range to the 
predicted instability ranges for the two pulsation models given 
in Lenz et al. (2008). Mode instability is mainly determined by 
the conditions in the stellar envelope; e.g., the high-frequency 
limit in S Scuti stars is very sensitive to the efficiency and treat- 
ment of envelope convection. In particular, near the terminal-age 
main sequence (TAMS), a small change in luminosity or temper- 
ature can imply very different stellar structures leading to differ- 
ent instability ranges. With a log g value of 3.6 + 0. 1 (Zima et al. 
2007), 44 Tau represents such an ambiguous case. 

The detection of a new frequency at 5.30 cd ' with an ampli- 
tude of 0.59 mmag poses a problem for the post-main sequence 
model. According to the results by Lenz et al. (2008), stability 
is predicted at such a frequency. We find that the low-frequency 
Hmit for instability cannot be shifted below this frequency by 
changing the input parameters (e.g., qtmlt) within reasonable 
limits. Consequently, there are two possibilities. 

(i) The main sequence model represents the correct model. 
This assumption requires a high metallicity of Z = 0.03 to fit the 



identified radial modes. However, spectroscopic measurements 
show solar abundances in the photosphere (Zima et al. 2007). 
Moreover, this model cannot reproduce the measured fundamen- 
tal parameters such as effective temperature. 

(ii) The post-main sequence scenario, which matches the 
fundamental parameters, is correct, but the computed envelope 
structure of the model does not describe the real physical con- 
ditions in the star perfectly. Therefore, the predicted range of 
unstable modes is narrower than the observed frequency range. 
This problem also exists for other 6 Scuti stars, such as FG Vir 
and for the two /? Cephei stars v Eri and 12 Lac, as stated by 
Dziembowski & Pamyatnykh (2008). These authors stress the 
necessity of an improvement in stellar opacity data by an opac- 
ity enhancement in the driving zone. Small opacity changes in 
the driving zone of 6 Scuti stars may also make the mode at 
5.30 cd"' unstable. 



7. Conclusions 

Two additional seasons of high-precision photometric data of the 
6 Scuti star 44 Tau were gathered with the Vienna Automatic 
Photoelectric Telescope in Arizona. A frequency analysis of the 
total data set including photometry from 2000 to 2006 led to 
detecting 49 pulsation frequencies, of which 15 are independent 
pulsation modes and 34 combination frequencies or harmonics. 

We find that the amplitudes of the combination frequencies, 
fj + fj, mirror the variations of the parent modes. The combi- 
nation parameter, fi, which relates the amplitudes of the combi- 
nation frequencies to those of the parent modes, is found to be 
different for different radial and nonradial parents. The combi- 
nation frequencies involving the radial fundamental mode have 
combination parameters, fi, of 0.003, while those involving the 
first radial overtone can be as high as w 0.009. 

The time-base of 5 years allows us to study amplitude vari- 
ability in more detail. Strong amplitude variability from year to 
year was found for the ( = I modes, while the two radial modes 
have essentially constant amplitudes. Moreover, the amplitudes 
of all of the ( = I modes decrease during most of the observed 
period. We examined several possible reasons for the amplitude 
variability of the ^ = 1 modes: 

(i) Precession of the pulsation axis: This causes variable am- 
plitudes because the mode visibility depends on the inclination 
angle. While the amplitudes of the radial modes remain constant, 
the amplitudes of /" = 1 modes will change. However, both ax- 
isymmetric and non-axisymmetric modes decrease at the same 
time; this observation cannot be explained in this scenario be- 
cause the visibility of m = modes behaves opposite to their 
non-axisymmetric counterparts. 

(ii) Beating of close frequencies: Due to the slow rotation of 
44 Tau, beating between the components of a rotationally split 
^ = 1 triplet has been examined within the limits of the measured 
rotation rate. Generally, the predicted beating timescales are too 
short compared to the observed amplitude modulation. However, 
beating with hitherto unresolved, additional close frequencies 
(which are not a component of the same rotational triplet) could 
produce amplitude modulation at the observed timescales. 

(iii) Resonance effects: The expected timescales for ampli- 
tude modulation in a 44 Tau model were computed. We find 
that only for a fraction of the variable { = I modes could the 
predicted variability timescales match the observed modulation 
periods. 

None of these three explanations can explain the observed 
amplitude variability completely. However, a combination of all 
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these effects cannot be excluded. Even more data would be re- 
quired to clarify this point. 

The newly found gravity mode at 5.30 cd"' extends the pre- 
viously known frequency range to lower frequencies. This detec- 
tion agrees better with the predicted instability ranges for main 
sequence models. However, our main sequence models cannot fit 
the measured fundamental parameters such as effective temper- 
ature. Our post-main sequence models provide a good fit of the 
fundamental parameters but predict stability at 5.30 cd"'. The 
problem with the present models underestimating the observed 
S Scuti star instability is common to many 6 Scuti stars and may 
be related to still incomplete opacity data. 
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